Magnetic dynamics with spin transfer torques near the Curie temperature 
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We use atomistic stochastic Landau-Lifshitz-Slonczewski simulations to study the interaction 
between large thermal fluctuations and spin transfer torques in the magnetic layers of spin valves. 
At temperatures near the Curie temperature Tc, spin currents measurably change the size of the 
magnetization (i.e. there is a longitudinal spin transfer effect). The change in magnetization of the 
free magnetic layer in a spin valve modifies the temperature dependence of the applied field-applied 
current phase diagram for temperatures near Tc- These atomistic simulations can be accurately 
described by a Landau-Lifshitz-Bloch + Slonczewski equation, which is a thermally averaged mean 
field theory. Both the simulation and the mean field theory show that a longitudinal spin transfer 
effect can be a substantial fraction of the magnetization close to Tc- 
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INTRODUCTION 



Spin transfer torque describes the interaction between 
the spin of itinerant, current-carrying electrons and the 
spins of the equilibrium electrons which comprise the 
magnetization of a fcrromagnet. This torque results 
from the spin-dependent exchange-correlation electron- 
electron interaction, and leads to the mutual precession 
of equilibrium and non-equilibrium spins around the total 
spin. In spin valves with sufficiently high current density, 
spin transfer torque can excite a free ferromagnetic layer 
to irreversibly switch between two stable configurations 
(typically along an easy-axis, parallel or anti-parallel to 
an applied magnetic field), or to undergo microwave os- 
cillations. Previous considerations of spin transfer torque 
mostly focus on the transverse response of the magneti- 
zation to spin currents [H [3j S 3 ■ This is appropriate 
since the temperatures used in spin valve experiments 
are substantially below the Curie temperature Tq of the 
ferromagnets, so that longitudinal fluctuations can be ig- 
nored. Near Tq, one expects an interplay between the 
large thermal fluctuations and the nonequilibrium spin 
transfer torque. Generally speaking, theories of critical 
phenomena in out-of-equilibrium systems have only re- 
cently been developed and there remain many open 
questions on this topic. 

Even far from the Curie temperature, temperature 
plays an important role in quantitatively analyzing the 
dependence of the magnetic orientation on the applied 
field and applied current. The effect of finite temperature 
on spin dynamics in the presence of spin transfer torque 
has been modeled the macrospin approximation (fixed 
magnetization length) by adding a Slonczewski torque to 
the Langevin equation describing the stochastic spin dy- 
namics [3, , and by solving the Fokker-Planck equation 
with the spin transfer torque term added to the deter- 
ministic dynamics [2]- The Keldysh formalism provides 
a formal derivation of the stochastic equation of motion 
for the non-equilibrium (i.e., current-carrying) sys- 
tem for a single spin of fixed magnitude. These treat- 
ments successfully describe the thermal characteristics 
of nanomagnets under the action of spin torques, such 



as dwell times and other details of thermally activated 
switching. 

For materials like GaMnAs, experiments are done near 
Tc , so that the size of the magnetization is substantially 
reduced from its zero temperature value (temperature 
in Kelvin), and undergoes sizeable fluctuations. In this 
case, the applicability of a macrospin model is not clear. 
For field-driven dynamics, there is theoretical work which 
accounts for longitudinal fluctuations near Tc This 
formal treatment culminates in the construction of the 
Landau-Lifshitz-Bloch equation (LLB), which is an ex- 
tension of the familiar Landau-Lifshitz equation with an 
additional longitudinal degree of freedom. In this work, 
we consider temperatures near the Curie temperature 
and include both longitudinal fluctuations of the mag- 
netization and the influence of spin transfer torque. 

There are a number of issues that complicate magnetic 
dynamics near Tc, including the temperature depen- 
dence of more basic magnetic properties such as magnetic 
damping and magneto-crystalline anisotropy, as well as 
the temperature dependence of the spin transfer torque 
itself. We use an atomistic approach for the stochastic 
dynamics of a local moment fcrromagnet with the inclu- 
sion of spin transfer torque. Such a model is more appro- 
priate for systems like the dilute magnetic semiconduc- 
tor GaMnAs. Our use of simple approximations for the 
temperature dependence of the magnetic anisotropy, de- 
magnetization field, and damping allow us to focus in the 
interplay between thermal fluctuations and spin transfer 
torque. We find that within this model, spin currents 
can change the size of the magnetization. We give an 
expression for this "spin-current longitudinal susceptibil- 
ity" , and propose an experimental scheme to measure 
this effect. 

We construct a Landau-Lifshitz-Bloch -I- Slonczewski 
(LLBS) equation to describe both longitudinal fluctua- 
tions and spin transfer torques. Following Ref. [T3 , we 
verify the applicability of the LLBS equation by compar- 
ing its results to the atomistic results. We then analyze 
the LLBS equation to find the applied field-applied cur- 
rent phase diagram for different temperatures. We find 
that critical switching currents are reduced by the same 
mechanism exploited in heat assisted magnetic recording. 
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namely the temperature-induced reduction in the mag- 
netic anisotropy [l^. We also find that regions of the 
phase diagram which have been experimentally unattain- 
able become relevant at high temperatures. The depen- 
dence of critical currents on temperature in these regions 
can provide quantitative details about the temperature 
dependence of spin transfer torque. 



II. METHOD 

To study the interplay between temperature and spin 
transfer torque, we consider a spin valve with a fixed layer 
magnetization in the -|-i-direction with Curie tempera- 
ture Tq, and a free layer with a smaller Curie tempera- 
ture Tq (see Fig. ([T])). This allows for a nearly temper- 
ature independent spin current flux incident on the free 
layer. We make the approximation that all of the incom- 
ing spin current is absorbed uniformly throughout the 
free layer magnetization. This approximation is based on 
two expectations. Substantial spatial and temporal in- 
homogeneities in the magnetization should induce rather 
irregular spatial patterns in the spin currents carried by 
propagating states. This will lead to large dephasing ef- 
fects, so that the total spin current should rapidly decay 
away from the interface as in the conventional picture 
of spin transfer torquesQ. In addition, in this tempera- 
ture regime, and for thin layers (« 3 nm) , magnetic non- 
uniformities in the direction transverse to current flow 
should be more substantial than non-uniformities along 
the current flow resulting from a localized spin transfer 
torque. 



I Fixed layer Free layer 




FIG. 1: Schematic of system, two ferromagnetic layers with 
different Curie temperatures. We suppose that Tq > Tq. 



A. Stochastic Landau-Lifshitz with spin transfer 

We adopt three approaches to model the system. The 
first is an atomistic lattice model of normalized spins S, 
which results in a stochastic Landau-Lifshitz equation 
(SLL). We include nearest-neighbor Heisenberg coupling 
with exchange constant J, and an easy-axis anisotropy 
field of magnitude TJan in the z-direction. To model the 
temperature dependence of the anisotropy, we make the 



ansatz that the magnitude of anisotropy at temperature 
T is proportional to the reduced magnetization m[T) = 
M,{T)/M^: 

H^,,{T) = H^,,{T = Q)m{T), (1) 

so that the anisotropy field on spin i is given by Hl^^ (T) = 
Han{T = 0)|S|S'f , where the bar indicates a spatial av- 
erage. A hard-axis anisotropy field with magnitude 
in the y-direction is added to model the demagnetization 
field of the thin layer. We make an ansatz for the form of 
this field to make the numerics more tractable. We take 
the demagnetization field to be uniform on all spins and 
given by i?^(T) = -HdiT = 0)Syy. This form of the 
hard- axis field ensures that ~ Ms{T), and roughly 
captures the non-local nature of the field. Finally, we 
include an applied field i/app iu the z-direction. The 
Hamiltonian for spin i is then: 



H, = J 



E 

jen.n. 



HdiT - 0)Sf (Sv 



(2) 



where the sum in the first term is over nearest neighbors, 
fiB is the Bohr magneton, and fio is the permeability 
of free space. To model nonzero temperatures, we add 
damping a and a stochastic field Ha to the equation of 
motion implied by Eq. with the standard statistical 
properties: 



(HSit)) = 0, 



a 2kBT 



(3) 

1 + .^ ,p 

where a, (3 are the Cartesian components of the field, ks 
is the Boltzmann constant, p is the magnetic moment 
on each lattice site, and 7 is the gyromagnetic ratio. 
We numerically integrate the equation of motion using 
a second-order Heun scheme . We add a Slonczewski- 
like spin transfer torque term to the equation of motion 
for the ith spin, which is given finally as: 

Si = -7Aio [Si X (Hcff + Hfi) - a (Si X Si X Hcfr) 

+Hi (S, X S, X z)] . (5) 

Hi parameterizes the spin transfer torque: Hj = 
—Ip^-b/ P'oe.^M^lA, where / is the applied current, p 
is the spin polarization of the current, is the zero 
temperature magnetization, £ is the free layer thickness, 
A is the transverse layer area, and — |e| is the electron 
charge. The effective magnetic field is given by Hcff = 
H^ppi + H.,.MS!z - (S^ y + J/ (mbA*o) E.en.n. S,- 
We use both a bulk geometry consisting of a iV = 48'^ 
periodic array of spins in 3 dimensions (simple cubic lat- 
tice), and a layer geometry with an array of 100 x 100 
X 15 spins. We employ the bulk geometry in compar- 
ing the stochastic model behavior with predictions from 
mean field theory, and the layer geometry for studying 
the effect of spin current on magnetization size. 
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B. Landau-Lifshitz-Bloch + Slonczewski equation 

In the second approach, we add a Slonczewski torque 
term to the LLB equation. To derive the LLB equa- 
tion, a probability distribution for the spin orientation 
is assumed, which is used to find the ensemble average 
of Eq. ([5]). In addition, the nearest neighbor exchange 
field is replaced by its mean-field value. The details of 
the derivation follow closely those in Ref. [iTI, so we omit 
them here. The final LLB+Slonczewski equation takes 
the form: 



m 



-7/^0 



2/l T 

(m X Hoff) -I- ^ _ ^ m • (aHoff -I- Hjz) m 



1 / k T 
— 7(1 ^ 1 m X m X (aHeff -|- Hiz) 



(6) 



with an effective field given by: 



H 



cff 



= H.. 



app-^ 



Ml 
2X 



1 m. 



(7) 



where is the zero temperature saturation magne- 
tization, m = M/Mg° is the dimensionless magnetiza- 
tion with magnitude between zero and one, me(T) is 
the zero field, zero current equilibrium magnetization: 
me(T) — B{jQ/kBT), and B is the Brillouin func- 
tion. x{T) is the longitudinal susceptibility: x{T) = 
{dme{T) /dHi^pp). Jq is the 0th component of the 
Fourier transformed exchange, and m is a vector with 
size between and 1. The spin transfer torque is param- 
eterized by Hj, as described in the previous section. The 
double cross product in Eq. ^ is the familiar Landau- 
Lifshitz damping term, which describes the relaxation of 
the magnetization direction to the nearest energy mini- 
mum. The term longitudinal to m distinguishes the LLB 
equation from the Landau-Lifshitz equation. This longi- 
tudinal term describes the relaxation of the size of the 
magnetization to its steady state value, which is deter- 
mined by the temperature, applied fields, and applied 
currents. 

The detailed dependence of the magnetic anisotropy 
on temperature is generally material specific. In our 
model, the anisotropy and demagnetization fields depend 
on temperature through their m dependence, and vary as 
ra^{T) and m{T), respectively. The magnetic exchange 
Jq can also depend on temperature. This dependence is 
stronger for ferromagnets with indirect exchange interac- 
tions (such as GaMnAs, where the magnetic interactions 
are mediated by hole carriers), and weaker for local mo- 
ment systems with direct exchange (such as Fe). For 
simplicity we treat Jq as temperature-independent. 

Finally we consider the standard Landau-Lifshitz equa- 
tion with a reduced but fixed saturation magnetization. 
We find in Sec. pil D|) that it is possible to appropriately 
modify the damping coefhcient in a standard Landau- 
Lifshitz approach so that the phase diagram it predicts 



agrees qualitatively with those predicted by the more 
complicated models. 



III. RESULTS 
A. Longitudinal spin current susceptibility 

In transition metal ferromagnets, longitudinal spin 
transfer, which is another way of saying spin accumu- 
lation, is typically quite small compared to the magne- 
tization and has a negligible effect on the magnetization 
dynamics. However, for temperatures close to the Curie 
temperatures, the longitudinal spin transfer can be a size- 
able fraction of the magnetization and can significantly 
affect the dynamics. 

Using the LLB+Slonczcwski equation, it is straightfor- 
ward to show that the change in the magnetization in the 
presence of spin current is 



(5m(/, T) = 



Hi x{T) 
MP 



a 



(8) 



This longitudinal spin transfer effect is demonstrated in 
Fig. [21 which shows the longitudinal susceptibility to 
magnetic field and spin current for a full stochastic sim- 
ulation with 100x100x15 spins. (In the figure, x is 
rescaled: the magnetic field is scaled by the exchange 
field Jq/^b/J-o, and the magnetization is scaled by M^.) 
In the simulation, the spins' polar angle is initialized to a 
uniform distribution between ^ = and 9 — ^max, where 
^max is chosen so that the initial spins' average is equal 
to the equilibrium value. We allow the system to relax to 
steady state, and find the value of the magnetization and 
its fluctuations by finding the average and standard de- 
viation over an interval of time (the appropriate time in- 
terval is temperature dependent). The fluctuations lead 
to the statistical uncertainty shown in Fig. 

The spin current susceptibility xi is defined as xi — 
Mj? (dm/dHj). We find that x and xio: correspond very 
well, demonstrating that Eq. ^ accurately describes the 
numerical stochastic model. 

The change in magnetization should be measurable. 
The fractional change in the magnetization compared to 
the zero-temperature saturation magnetization is 



Sm — 



'X(T) 



(9) 



For T/Tc = 0.95, so that {x ■ Jo/ fiofisMs) = 7 (from 
Fig. (12)), and with an exchange field of Jo//^b/^o — 
1.2 X 10^ A/m (which corresponds to a Tc of 150 K 
in a cubic nearest neighbor Heisenberg model), = 
10^ A/m, I/A = 10"A/m2, p = 0.5, a = 0.01, and 
£ ^ 3 nm gives a change compared to the zero tem- 
perature value of Sm =5 %. Since the magnetization 
is reduced to approximately 20 % of its zero tempera- 
ture value at T/Tc = 0.95, the fractional change in the 
magnetization is approximately 25 %. 
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FIG. 2: The magnetic field and spin current susceptibility 
versus temperature for the stochastic Landau-Lifshitz equa- 
tion in the layer geometry. The spin current susceptibility 
is multiphed by a. The error bars indicate statistical uncer- 
tainty (one standard deviation). In the plot, x is rescaled by 
/ioPsM°/J. 



A notable aspect of this longitudinal spin transfer is 
that the size of the magnetization can either be increased 
or decreased according to the direction of current flow. 
For electron flow from fixed to free layer, the free layer 
moment increases, while electron flow in the opposite 
direction decreases the free layer moment. This con- 
trasts with current-induced Joule heating, which always 
decreases the magnetization. 

This distinction can be exploited to probe the longi- 
tudinal spin transfer by using the experimental scheme 
shown in Fig. ([3]). We consider the case where ^ T > 
Tq. We choose sign conventions such that a positive i?app 
aligns with the fixed layer, and a positive current repre- 
sents electron flow from fixed to free layer. In the absence 
of a longitudinal spin transfer (x/=0, black line in Fig. 
([3])), the application of a magnetic field will partially or- 
der the free layer to align or anti-align with the fixed 
layer. This should cause the resistance R of the device to 
change in some way, according to the giant magnetoresis- 
tance effect and magnetic order induced in the free layer 
(the detailed dependence of R on i?app is not important 
here). If a positive current is applied, then the longi- 
tudinal spin transfer induces partial ordering of the free 
layer, so that m (iJapp = 0) = +xiHjo/M°. Then the 
curve of TO(iJapp), and therefore the curve R (^/^app) is 
simply shifted by +XiHjo /x (the red dashed curve in 
Fig. ([3])). If a negative current density — is applied, 
then m (i?app = 0) = —xiUjo /M^ and the m (J^app) and 
R{Ha,pp) curves are shifted by —XiHio/x (black dotted 
curve in Fig. ([3])). This shift represents a unique signa- 
ture of longitudinal spin transfer. 

Using the same parameters as before, we estimate a to- 
tal shift 5 = 2xiHjo/x between R (-ffapp) for positive and 



negative current to be on the order of 8 x 10^ A/m (« 
1 T). Eq. ^ indicates that materials with small ex- 
change field (or small Ic), and those that can support 
large current densities show the effect most strongly. 
This suggests that weak metallic ferromagnets such as 
Gd (Tc = 300 K), and Fe alloys such as FeS2 and 
FeBes (Tc = 270 K) [ll| may be good candidates for 
free layer material. 
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FIG. 3: Experimental scheme for detecting longiudinal spin 
transfer: for S> T > T^, an apphed field Happ changes the 
resistance R via the magnetoresistance efi'ect. The applica- 
tion of a positive and negative current density of magnitude 
/" shifts m(i/app) in the positive and negative direction, re- 
spectively, via longitudinal spin transfer. The -R(iifapp) curves 
therefore shift to the positive and negative directions. 



B. Landau-Lifshitz-Bloch-Slonczewski vs 
Stochastic Landau-Lifshitz 

In this section, we compare the results obtained from 
the full 3-dimcnsional stochastic LL+S equation with 
those obtained from the mean-field LLBS equation. In 
our numerics, we rescale time t to as r = (7J//iB)i, 
which rescales the magnetic fields _ffoff by the exchange 
field 7?ox = J/^J'Q^J'B■ Dimensionless fields are denoted by 
lowercase: /lapp = ^^appMoMB/^/, etc. The dimensionless 
spin torque is denoted by japp, where japp = HjfiQ^B/J- 
We consider a current-induced magnetic excitation for 
the bulk lattice geometry at various temperatures. The 
average magnetization is initialized at a 45° angle with 
respect to the -l-z-direction (the individual spins' initial 
direction is distributed uniformly within 3° in the 9, 4> 
direction about = 45°, (/) = 0°). The spin transfer 
torque is applied to excite the magnetization away from 
the z-direction. The parameters used are an applied field 
of /lapp = 0.0001, a demagnetization field of /id = 0.01, a 
current of Japp — 

-0.0002, and damping of a = 0.1 (the 
artificially high damping was chosen to allow the numer- 
ical simulation of the switching to be carried out in a 
reasonable time). The time step used for the numerical 
integration is dr — 0.0002. We vary the temperature T, 
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and present results in terms of the scaled temperature 
T = T/Tq. 

As we increase temperature, we obtain trajectories of 
varying complexity. Fig. ^ compares the LLBS and sev- 
eral realizations of the stochastic Landau-Lifshitz equa- 
tion. For this range of parameters, the magnetic dynam- 
ics evolves from steady oscillations to current induced 
switching as the temperature is increased. Generally, the 
level of correspondence between the two is qualitatively 
good, although it varies between different realizations of 
the stochastic dynamics. We can conclude from this data 
that the LLBS equation qualitatively captures the fea- 
tures of the full stochastic simulations. 

The trajectories for t = 0.08 indicate that a realization 
of stochastic dynamics can exhibit the crossover from pre- 
cession to stable switching, whereas at this temperature 
the trajectory obtained with the LLBS equation shows 
only oscillations. This illustrates an important distinc- 
tion between the stochastic Landau-Lifshitz and LLBS 
models. The LLBS is an equation for the thermally av- 
eraged magnetization, derived using an assumed proba- 
bility distribution function (in this case, a distribution 
function most appropriate for temperatures well above 
and below energy barriers). For this reason, the LLBS 
does not contain information about fluctuations, and in 
particular does not capture stochastic switching over the 
energy barrier. The fluctuations may be obtained by 
solving the Fokker- Planck equation, or by supplement- 
ing the LLBS equation with stochastic fields, as done in 
Ref. [H 

C. Applied field-applied current phase diagram 

Both high temperatures and the longitudinal degree of 
freedom change the applied field-applied current phase 
diagram of the free magnetic layer. Fig. ([5]) shows the 
generic topology for regions of stability for the parallel 
("P", or -t-z-direction) and antiparallel ("AP", or —z- 
direction) fixed points. We focus on the stability of the 
AP configuration for positive applied fields (the dashed 
boundary in the upper-half- plane of Fig. We first 

briefly describe the main qualitative features before pro- 
viding a mathematical description. For applied fields be- 
tween hanm^ and h^nm^ + ham, the stability boundary 
is a horizontal parabola, while for other values of ap- 
plied field, the stability boundary is linear with slope 
1/a. For applied fields with magnitude less than h^n^, 
there is hysteresis in the current switching. For T — 0, 



this phase diagram reduces to the known form found ex- 
perimentally As T increases, the size of the hys- 
teretic region (and the switching current) decreases. Also 
the range of field with the parabolic boundary decreases, 
and the outer edge of the parabola gets pulled in closer 
to 0. For sufficiently high temperatures, this parabolic 
stability boundary should be experimentally accessible. 

A quantitative description of the phase diagram follows 
from Eq. © . We determine the stability of fixed points 
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FIG. 4: Comparison of the (i-component) magnetization time 
evolution with spin transfer torque for the atomistic stochas- 
tic simulation and the LLB+Slonczewski equation for various 
reduced temperatures T' — T/Tc- The dashed line gives the 
LLB-|-Slonczweski trajectory, while the solid lines show vari- 
ous realizations of the stochastic trajectory. The dimension- 
less time T is given by r = {-yj/ fis) t. 



using the standard method of linearizing Eq. ^ about 
a fixed point and finding parameter-dependent eigenval- 
ues A. A positive real part of A indicates a loss of sta- 
bility. This analysis leads to the following condition 
for instability of the antiparallel configuration (where 
it should be noted that m depends on japp through 

TO — rUe + X (ji + , and x is the rescaled suscep- 
tibility, given by X = X (Jq/ f^ofJ-BM°)): 
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This leads to a cubic equation for j™*. Assuming nie > x(^+— )- and expanding to 0th order in x 
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leads to an approximate, closed form for j'app- Again we 



distinguish between different regimes of applied field. For 



h 1 

li H nip 

2 



; 1 - ST' 
1 -T' 



(10) 



where again m,, is the equilibrium magnetization in the 
absence of applied field and applied current. Eq. JTU]) 
shows that the slope of the boundary is temperature in- 
dependent, and is given by 1/a (the intrinsic damping a 
is assumed to be temperature independent). The tem- 
perature independence of the slope follows from the fact 
that the spin transfer torque increases like l/m(T), but 
the effective damping rate also increases as l/m(T). The 
intercepts of this boundary line are temperature depen- 
dent due to the temperature dependence of m. The con- 
tribution from the easy-axis anisotropy field has an ad- 
ditional temperature dependence, but the magnitude of 
this field is much smaller than the demagnetization field, 
so it does not play an important role. The critical current 
at zero field is reduced by m(T) because of the reduction 
in the demagnetization field. This is important because 
the demagnetization field is usually larger than applied 
fields, and is therefore the primary impediment to cur- 
rent induced switching. Its reduction through increased 
temperature offers a route to reduced critical switching 
currents. 



For /la 



< h < 



hdin, a very large spin 



torque is required to stabilize the AP configuration. The 
values of current for which the AP configuration is stabi- 
lized are much higher than those attainable experimen- 
tally, so that for this range of fields the AP configuration 
is not seen [18]. The approximate critical current along 
the AP stability boundary is: 



'.\/h{h, 



drrie 



h) 



./app 



1 -T' 



(11) 



The reduction in the outer boundary of the parabolic 
stability line is reduced at high temperature, and this re- 
duction can also be traced back to the reduced magnetic 
anisotropy. For low temperatures, the application of spin 
transfer torques results in a elliptical precession mostly in 
the easy plane about the —z fixed point. To stabilize the 
AP configuration in this regime, the spin transfer torque 
must overcome the precessional torque (usually, the spin 
transfer torque must overcome the much weaker damp- 
ing torque). Assuming h — h(im/2 for definiteness, the 
precessional torque decreases with T as hdm(T), while 
the spin transfer torque increases like 1/m. This im- 
plies a value for the maximum reach of the parabola of 
japp = m2(r)/id/(2(l - T)). Plugging in typical value s 
for material parameters (the same used in Sec. (jlll Ap ) 
leads to a critical current of IQi^A/m^ for T = 0.95Tc. 
This is an order of magnitude smaller than the zero tem- 
perature case. The behavior of this critical current versus 
temperature at a fixed applied field is shown in Fig. 
(Solid line gives LLBS result). It should also be noted 



that the stochastic trajectories (shown in Fig. (|4])) in- 
dicate that thermal fluctuations can effectively drive the 
system out of the precessional state and into the static 
antiparallel configuration. 




FIG. 5: Schematic of parallel/anti-parallel stability versus ap- 
plied field and applied current. The hysteretic box near the 
origin and the fully unstable regions (white parabolic shapes) 
contract in size with increasing temperature. 



D. Comparison with Landau-Lifshitz-Slonczewski 

The Landau-Lifshitz-Slonczewski (LLS) equation can 
be modified to emulate the LLBS equation. Based on 
the qualitative behavior of the LLBS equation, a suit- 
able form for a temperature dependent LLS equation for 
a nanomagnet of reduced magnetization size m and ori- 
entation h is; 



7/io fi X Hcff n X n X H, 



Hi 



eff 



-n X n X z 



where Hp 



loff = Happ - rnHdUyy + m-^Hi,^nzZ, and the 
temperature dependence is contained entirely in m(T). 
Clearly the divergence of the damping at T = Tc is un- 
physical, however a more detailed treatment of damping 
near Tq is beyond the scope of this paper. The differences 
between this LLS equation and the LLBS equation are 
quantitative (as opposed to qualitative) in nature. One 
difference is in the dependence of the critical current on 



temperature for 



< h < /lanW^ 



hd-m. Fig. dSl) 



shows the prediction based on the LLS equation. 

The LLS equation equation neglects the longitudinal 
spin transfer and applied field susceptibility, which are 
responsible for dynamically changing the size of the mag- 
netization (and therefore the size of the effect fields) dur- 
ing a switching event, or other magnetization dynamics. 
However, Fig. ^ shows qualitative agreement between 
the critical currents found in both LLBS and LLS mod- 
els. This is indicative of the fact that for the applied 



7 









LLBS 

---LLS 


X 

X 


X 

X \ 
X \ 
X \ 
X \ 
X \ 

X \ 

\ \ 
\\ 
\\ 



0.2 0.4 0.6 0.8 1 
T/T 

c 



FIG. 6; Critical current versus temperature for LLBS and 
LLS equations. The parameters are: /lapp ~ —0.001, /id ~ 
0.01, /lan ~ 0.0001. Recall that all fields are scaled by the 
exchange field. 



field-applied current phase diagram, the spin-current and 
applied-field longitudinal susceptibilities play a role that 
is secondary to the more pronounced effects of tempera- 
ture reduced anisotropics. 

IV. DISCUSSION 

Spin transfer torques can affect the longitudinal fluctu- 
ations of a ferromagnet near its critical temperature. To 
consider these effects, we studied an atomistic, stochastic 
Landau-Lifshitz-Slonczewski simulation at high temper- 
atures. We find that there is a longitudinal spin trans- 
fer effect, and estimate that at temperatures near Tq, 
spin currents can measurably change the size of the mag- 
netization. We then supplemented the Landau-Lifshitz- 
Bloch equation with a Slonczewski torque term, and ver- 
ified that this model captures the qualitative features 
of the stochastic simulations. We showed that the ap- 
plied field-applied current phase diagram undergoes large 
changes in the presence of high temperatures, and that 
these changes may be useful for reducing critical switch- 
ing currents and for studying the detailed behavior of the 



temperature dependence of the spin transfer torque. It 
should be emphasized that these results are predicated 
on a disordered local moment model of a ferromagnetic 
phase transition. This model leads to an effective damp- 
ing that increases with temperature as l/m{T), which ef- 
fectively counteracts the similar l/m(T) increase in the 
magnitude of spin transfer torque. Materials that un- 
dergo a Stoner transition should also have a l/m{T) de- 
pendence for the spin transfer torque, but a different tem- 
perature dependence for damping. Such materials should 
therefore behave differently than the model considered 
here. 

The experimental system relevant for the effects we 
describe (shown schematically in Fig. ([1])) should be rel- 
atively straightforward to fabricate. Jiang et al. con- 
sidered a similar system [19j, although that work dealt 
with other issues such as the ferrimagnet compensation 
point for magnetization and total angular momentum. 
By considering simpler ferromagnets with different Curie 
temperatures, the role of temperature may be more eas- 
ily inferred. It is of course necessary to account for Joule 
heating in assessing the detailed temperature dependence 
of the spin transfer torque. However recent experiments 
on domain wall motion illustrates the feasibility of com- 
pensating for this effect On the other hand, ex- 
periments conducted at fixed current with varying am- 
bient temperatures and applied fields may offer a more 
straightforward route to observing the longitudinal spin 
transfer effect. 

Many experiments done with dilute magnetic semi- 
conductors deal with domain wall motion, where ther- 
mal effects play an important role in even the quali- 
tative aspects of the domain wall behavior [Soj. There 
are additional challenges associated with extending this 
work from spin valves to continuous magnetic textures. 
Among these is the renormalization of the exchange in- 
teraction associated with the coarse graining of the mag- 
netization, which becomes more important at higher tem- 
peratures [2l| . In addition, the crucial role played by the 
demagnetization field in intrinsic domain wall pinning 
implies that the finite temperature treatment of the de- 
magnetization field must also be handled more carefully. 
For these reasons the spin valve geometry may provide 
greater experimental control and admit a simpler theo- 
retical description. 
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